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Spacetime foam can be modeled in terms of nonlocal effective interactions in a classical nonfluctu- 
ating background. Then, the density matrix for the low-energy fields evolves, in the weak-coupling 
approximation, according to a master equation that contains a diffusion term. Furthermore, it is 
argued that spacetime foam behaves as a quantum thermal field that, apart from inducing loss 
of coherence, gives rise to effects such as gravitational Lamb and Stark shifts as well as quantum 
damping in the evolution of the low-energy observables. These effects can be, at least in principle, 
experimentally tested. 
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I. INTRODUCTION 

It seems natural to assume that spacetime at the 
Planck scale must have a very complicated and ever- 
changing topology. Indeed, it was Wheeler [Q who first 
suggested the foamlike structure of spacetime as 
an inescapable ingredient of the yet-to-be-built quantum 
theory of gravity. Since then, various spacetime foam 
components have been proposed: wormholes jJJ^, vir- 
tual black holes ||^, and quantum time machines |^,^ 
among them. 

The quantum theory of gravity suffers from problems 
that have remained unsolved for many years. They 
originate in the fact that gravity deals with the frame 
in which everything takes place, i.e., with spacetime, 
in sharp contrast with any other interaction, for which 
spacetime is a passive frame. When gravity is brought 
onto the scene, the frame itself becomes dynamical. It 
suffers the quantum fluctuations of the other interactions 
and, even more, introduces its own fluctuations, thus be- 
coming an active agent in the theory. 

We are used to putting everything into spacetime, so 
that we can name and handle events. General relativ- 
ity made spacetime alive and in this sense, was a major 
change. But, although dynamical, the relations between 
different events were still sharply defined. Quantum me- 
chanics changed this, too. In such a dynamical frame, 
objects became fuzzy; exact locations were substituted 
by probability amplitudes of finding an object in a given 
region of space at a given instant of time. 

A quantum uncertainty in the position of a particle 
implies an uncertainty in its momentum and, therefore, 
due to the gravity-energy universal interaction, would 
also imply an uncertainty in the geometry, which in turn 
would introduce an additional uncertainty in position of 
the particle. The geometry would thus be subject to 
quantum fluctuations that would constitute the space- 
time foam and that should be of the same order as the 
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geometry itself at the Planck scale. This would give rise 
to a minimum length [ pO[ beyond which the geometri- 
cal properties of spacetime would be lost, while on larger 
scales it would look smooth and with a well-defined met- 
ric structure. 

The quantum structure of spacetime would be relevant 
at energies close to Planck scale and one could expect 
that the quantum gravitational virtual processes that 
constitute the spacetime foam could not be described 
without knowing the details of the theory of quantum 
gravity. However, the gravitational nature of spacetime 
fluctuations provides a mechanism for studying the ef- 
fects of these virtual processes in low-energy physics. In- 
deed, virtual gravitational collapse and topology change 
would forbid a proper definition of time at the Planck 
scale. More explicitly, in the presence of horizons, closed 
timelike curves, topology changes, etc., any Hamiltonian 
vector field that represents time evolution outside the 
fluctuation would vanish at points inside the fluctuation. 
This means that it would not be possible to describe the 
evolution by means of a Hamiltonian unitary flow from 
an initial to a final state and, consequently, quantum 
coherence would be lost. These effects and their order 
of magnitude would not depend on the detailed struc- 
ture of the fluctuations but rather on their existence and 
global properties. In general, the regions in which the 
asymptotically timelike Hamiltonian vector fields vanish 
are associated with infinite redshift surfaces and, con- 
sequently, these small spacetime regions would behave 
as magnifiers of Planck length scales transforming them 
into low-energy modes as seen from outside the fluctu- 
ations Therefore, spacetime foam and the related 
minimum length would affect low-energy physics, so that 
low-energy experiments would effectively suffer a non- 
vanishing uncertainty. In this situation, a loss of quan- 
tum coherence would be almost unavoidable jlj]. In fact. 
Hawking ^ has pointed out that scalar fields may lose 
coherence extremely fast and that the loss of quantum 
coherence might also be responsible for the vanishing of 
the 9 angle of quantum chromodynamics. 

In this paper, we show that spacetime foam behaves as 
a quantum thermal bath with a nearly Planckian temper- 
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aturc that has a weak interaction with low-energy fields. 
As a consequence, other effects, apart from a loss of co- 
herence, such as Lamb and Stark transition-frequency 
shifts, quantum damping, and cold diffusion, character- 
istic of systems in a quantum environment nat- 
urally appear as low-energy predictions of this model. 
A brief account of these results has already appeared in 
Ref . . This kind of quantum gravitational effects can 
be, in principle, experimentally tested (see, e.g., Refs. 
15 1^), as we also argue in this work. 

This paper is organized as follows. In Sec. ||, we pro- 
pose an effective model of spacetime foam in terms of 
nonlocal interactions and argue that it is equivalent to 
a local theory with a stochastic classical Gaussian noise 
source. In Sec. [II , a master equation for the low-energy 
density matrix is obtained and the diffusion term is stud- 
ied. Section [V is devoted to the quantum effects that 
have not been taken into account in previous sections and 
derive a master equation that includes them. It is also 
shown that spacetime foam can be described as a quan- 
tum thermal bath and the consequences of this effective 
behavior are analyzed. We close this section with a short 
discussion about the kind of experiments and observa- 
tions that could be sensitive enough to test these effects. 
In Sec. ^ we study the role that some of the compo- 
nents of spacetime foam (wormholes, virtual black holes 
and quantum time machines) play in the effective theory. 



We summarize and conclude in Sec. VI 



II. EFFECTIVE INTERACTIONS 

In this section, we will construct an effective theory 
for the evolution of low-energy fields in spacetime foam, 
where we possibly have a finite resolution limit because 
the notion of distance is not valid at the quantum gravi- 
tational scale. 

With this aim, we will substitute spacetime foam by 
a fixed classical nonfluctuating background with low- 
energy fields living on it. We will perform a 3+1 foli- 
ation of the effective spacetime that, for simplicity, will 
be regarded as fiat, t denoting the time parameter and 
X the spatial coordinates. Spacetime foam features, i.e., 
the gravitational fluctuations and the minimum length 
generated by them, will be characterized by nonlocal in- 
teractions. They will relate spacetime points that are suf- 
ficiently close in the effective nonfiuctuating background, 
where a well-defined notion of distance exists. These ef- 
fective nonlocal interactions will be described in terms of 
local interactions as follows. 

Let us consider a basis {hi{t)} of local gauge- invariant 
interactions at the spacetime point {x,t), each element 

consisting of factors of the form Z^"'^^*'"** ['/'(a:, ^)]^", and 
(j) being the low-energy field strength of spin s. As a 
notational convention, each index i implies a dependence 
on the spatial position x by default; whenever the index 
i does not carry an implicit spatial dependence, it will 



appear underlined i. Also, any contraction of indices 
(except for underlined ones) will entail an integral over 
spatial positions. Then, the nonlocal effective interaction 
can be included in the Euclidean action by means of a 
term of the form 



N 



where Lnt is the A^-local interaction term. 



(2.1) 



Ijv = ^ j dh---dtNc''""''{ti...tN)h,,{h)---h,^{tN). 

(2.2) 

The dimensionless functions c*^ ' *"(ti . . .ijv) cannot de- 
pend on the location of the gravitational fluctuation itself 
because of conservation of energy and momentum: the 
fluctuations do not carry energy, momentum, or gauge 
charges. Thus, diffeomorphism invariance is preserved, at 
least at low-energy scales, provided that the coefficients 
c*i "*"(ti . . . Im) only depend on relative positions. This 
invariance cannot be expected to hold at the Planck scale 
as well. However, this violation of energy-momentum 
conservation is safely kept within Planck scale limits |l^ , 
where the processes will no longer be Markovian. 

The coefficients c*^ ' '"(ii . . .t^) must vanish for rel- 
ative spacetime distances larger than the length scale r 
of the gravitational fluctuations. Indeed, if the gravi- 
tational fluctuations are smooth in the sense that they 
only involve trivial topologies or contain no horizons, the 
coefficients c*i ' '"(ii . . Am) will be iV-point propagators 
which, as such, will have infinitely long tails and the 
size of the gravitational fluctuations will be effectively 
inflnite. In other words, we would be dealing with a 
local theory written in a nonstandard way. The gravi- 
tational origin of these fluctuations eliminate these long 
tails because of the presence of gravitational collapse and 
topology change. This means that, for instance, virtual 
black holes will appear and disappear and horizons 
will be present throughout. As Padmanabhan ||l|] has 
also argued, horizons induce nonlocal interactions of fi- 
nite range since the Planckian degrees of freedom will be 
magnified by the horizon (because of an infinite redshift 
factor) thus giving rise to low-energy interactions as seen 
from outside the gravitational fluctuation. Virtual black 
holes represent a kind of components of spacetime foam 
that because of the horizons and their nontrivial topol- 
ogy will induce nonlocal interactions but, most probably, 
other fluctuations with complicated topology will warp 
spacetime in a similar way and the same magnification 
process will also take place. 

Finally, the coefficients c'^ ' *" [ti . . . t^) will contain a 
factor [e"'^*^''^/^]^, S{r) being the Euclidean action of the 
gravitational fiuctuation, which is of the order (r/Z*)^. 
This is just an expression of the idea that inside large 
fiuctuations, interactions that involve a large number of 
spacetime points are strongly suppressed. As the size of 
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the fluctuation decreases, the probabihty for events in 
which three or more spacetime points are correlated in- 
creases, in close analogy with the kinetic theory of gases: 
the higher the density of molecules in the gas, the more 
probable is that a large number of molecules collide at 
the same point. The expansion parameter in this exam- 
ple is typically the density of molecules. In our case, the 
natural expansion parameter is the transition amplitude. 
It is given by the square root of the two-point transition 
probability which in the semiclassical approximation is 
of the form e"'^'^''^. 

A simple calculation shows that 

Tn - 6^(r/0^^-^ {{{hllf-^-^'^^^^-^ (2.3) 

1=1 

where e = e~^^^^^'^{r/U)'^ . Indeed, In contains a factor 
|^g-S(r)/2jAr QQjjjjjj^g from the coefficient c*i ' '"(ti . . .t^r), 
as discussed above; each interaction hi provides a factor 
there are also N integrals over space- 
time positions, iV — 1 of which are integrals over relative 
positions and therefore give a factor r'^ each; and, finally, 
the integral over the global spacetime position provides 
an additional factor l^. The interaction term Tjy has 
contributions from three different length scales, Planck 
length 1*, the size of the gravitational fluctuations r, and 
the low-energy length scale through the ratios r/U, 
r/l and l/h: the factor e^(r/Z)^^~** depends only on 
the first two and is common to all powers n and spins 
s while the factors (Z*/Z)^"i(^+''i^~^ depend only on the 
low-energy scale (in Planck units) and contain the infor- 
mation about the different kind of interactions involved. 

The contributions of the trilocal and higher effective 
interactions are, at most, of order e'^. Therefore, in the 
weak-coupling approximation, i.e., up to second order in 
the expansion parameter e, they can be ignored. On the 
other hand, the local terms Xq and Xi can be absorbed 
in the bare action. Indeed, the coefficient c appearing in 
Xq is constant; the coefficients d{t) in Xi cannot depend 
on spacetime positions because of diffeomorphism invari- 
ance and are therefore constant as well. Consequently, we 
can write the nonlocal interaction term in the Euclidean 
action as the bilocal contribution 

lint - I Jdtdt'c'^t - t')h,{t)h,{t'), (2.4) 

where we have renamed c*-'(i, t') as &^{t — t'). This coeffi- 
cient is symmetric in the pair of indices ij and depends on 
the spatial positions Xi and Xj only through the relative 

distance jxi — x^ j. It is of order e"'^^''' and is concentrated 
within a spacetime region of size r. 

The effect of a single spacetime fiuctuation can be de- 
scribed in the path integral approach by adding a con- 
tribution J'D(j)e~^° lint to the bare low-energy Euclidean 
path integral / 'D(j>e~^° , Iq being the bare low-energy ac- 
tion. If we consider J\f indistinguishable gravitational 



fiuctuations, the contribution is JVcjie"^" (Lmt)''^ /J^^-- 
Thus, summing over any number M of them, we obtain 
the path integral /X>0e"-^°+^'°' . 

The bilocal effective action above does not lead to a 
unitary evolution for the low-energy fields because there 
exist different trajectories that arrive at a given conflgu- 
ration (0, 0); the future evolution depends on these past 
trajectories and not only on the values of (f> and (/) at that 
instant of time. Therefore, it is not sufficient to know 
the fields and their time derivatives at an instant of time 
in order to know their values at a later time: we need to 
know the history of the system, at least for a time r. As 
a consequence, the system cannot possess a well-defined 
Hamiltonian vector field and undergoes an intrinsic loss 
of predictability . 

The exponential of the interaction term e^'"' can be 
written as [ pO| 

jvae'^ Jdtdt''rijit-t')a'it)a^t')^- Jdta'{t)hi{t) ^ ^2.5) 

where, the continuous matrix "fij(t — t') is the inverse of 
c'J(t-t'), i.e., 

j dt"j,k{t - <")c"^(<" - t') Si6(t - t'). (2.6) 

Note that the quadratic character of the distribution for 
the fields a* is a consequence of the weak-coupling ap- 
proximation (second order in e), which keeps only the 
bilocal term in the action. Beyond the weak-coupling 
approximation, higher-order terms would introduce de- 
viations from this noise distribution. Note also that we 
have a different field a* for each kind of interaction hi. 
Thus, we have transferred the nonlocality of the low- 
energy fields (j> to the set of fields a*, which arc non- 
trivially coupled to it. 

If we now perform a Wick rotation back to Lorentzian 
spacetime, we see that the path integral has the form 

JvaP[a] J-p^e'^So+Jdta^WhAt)]^ ^2.7) 

where Sq is the low-energy Lorentzian action and 
P[a] ^ / dtdt''ri,it-t')a' {t)a^t') ^- J dta\t)hi{t) ^2,8) 

is the Gaussian probability distribution with correlation 
functions c^^ {t — t') for the stochastic nonlocal fields a* 
that represent spacetime foam and which are not affected 
by the Wick rotation. 



III. CLASSICAL DIFFUSION 

The analysis of the previous section ignores in a way 
the quantum nature of gravitational fluctuations such as 
virtual black holes or quantum time machines. Indeed, 
the fields a* represent quantum gravitational spacetime 
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foam but, as we have seen, the path integral for the whole 
system does not contain any trace of the dynamical char- 
acter of the fields a*. It just contains a Gaussian prob- 
ability distribution for them. The path integral above 
can then be interpreted as a Gaussian average over the 
classical noise sources a*. Classicality here means that 
we can keep the sources a* fixed, ignoring the noise com- 
mutation relations in a kind of zeroth-order semiclassical 
approximation, and, at the end of the calculations, we 
just average over them. The next section will be devoted 
to the quantum noise effects generated by spacetime foam 
that we are ignoring here. 

The master equation that governs the Lorentzian dy- 
namics of the low-energy fields in foamlike spacetimes is 
derived in what follows. For each fixed set of fields a*, the 
evolution equation for the density matrix pa (t) , obtained 
with the Hamiltonian 

H^{t) = Ho + a\t)h,, (3.1) 

Hq being the bare Hamiltonian of the low-energy field, is 

Po.{t) ^ ~i[Ho, Pa{t)] - ia\t)[h,, pa{t)]. (3.2) 

In the interaction picture, this equation becomes 

p\{t) = -za\t)[h\{t),p\{t)i (3.3) 

where 

pLW = U+{t)pa{t)Uo{t), (3.4) 

h\{t) = U+it)h,Uo{t) (3.5) 

with Uo{t) = e"*^"*. Averaging this equation over a* 
would provide a master equation for the density matrix 
of the system although it would not be very useful since 
it would contain terms of all orders in a*. To avoid this 
problem, we integrate this equation between an initial 
time to and t: 

Pait) = P'M - ^ fdt'a\t')[h\{t'),p'^{t')] (3.6) 

Jta 

and introduce this formal solution back to the differ- 
ential equation for p\ noting that at the initial time 
pL(^o) — P^{to) does not depend on a* (if this were not 
the case, there would be an extra renormalization term 
in the Hamiltonian): 

p'^{t) = -ta\t)[h\{t),p\to)] 

- fdt'a\t)a\t')m),[h){t'),p\{t')]]. (3.7) 

J to 

Next, we perform the Gaussian average over a* taking 
into account that ^^(t) does not depend on a* at zeroth 
order but only at first order, i.e., p\{t) = p^{t) + 0{a) 
with p^{t) = {Pa{t)) and keep terms up to second order 
in e (weak-coupling approximation). We then obtain the 
following equation for p^{t): 



p\t) = - f dT{a\t)a^{t^T)) 

JO 

x[h\{t),[h){t^T),p\t~T)]], (3.8) 

where we have made a change of integration variables 
from t' to T = t — t' . We also assume that p^{t) hardly 
changes within a correlation time r (Markov approxi- 
mation), so that p^{t — r) ~ p'(0- This amounts to 
ignore terms of order in the master equation. The 
initial condition can be taken at io = — oo, so that 
the integration range is now (0, oo). Note that, since 
{a^{t)a^{t — t)) = d-'ir) is nonvanishing only for t < r, 
this limit to —oo just implies that the evolution must 
take place over periods of time much larger than the cor- 
relation time r for the approximation to be valid. The 
resulting master equation in the interaction picture is 
then 

poo 

p\t) = / dTc}={T)[h\{t), [h){t-T),p\t)]l (3.9) 
JO 

Transforming this equation back to the Schrodinger 
picture, we obtain the equation 

/•oo 

p = -i[Ho,p\- drc'^T) [h,, [h){-T),p]] . (3.10) 
Jo 

Since h){-T) = {-T)hjUo{-T) and Uo{t) = 1 + 
0{t/1), the final form of the master equation for a 
low-energy system subject to gravitational fluctuations 
treated as a classical environment and at zeroth order in 
r/l (the effect of higher order terms in r/l will be thor- 
oughly studied in the next section) is pl| 

POO 

p^-i[Ho.p]~ drc*^(r)[/i„[/i,,p]]. (3.11) 
Jo 

The first term would also be present in the absence of 
fluctuations, since it governs the low-energy Hamiltonian 
evolution. The second term is a direct consequence of 
the foamlike structure of spacetime and the related exis- 
tence of a minimum length. It is a diffusion term which 
will be responsible for the loss of coherence. Note that a 
dissipation term, necessary to preserve the commutation 
relations under time evolution, is not present. However, 
we have considered the classical noise limit, i.e., the fields 
a* have been considered as classical sources and the com- 
mutation relations are automatically preserved. We will 
see that the dissipation term, apart form being of quan- 
tum origin, is r/l times smaller than the diffusion term 
and we have only considered the zeroth order approxi- 
mation in r/l. 

The diffusion term induces a characteristic decoherence 
time Td that can be easily calculated. Indeed, the interac- 
tion Hamiltonian density hi is of order /~''(Z*/Z)^"i^^"'"''i^ 
and d'-' {t) is of order e"'^^'"-'. Furthermore, the diffusion 
term contains one integral over time and two integrals 
over spatial positions. The integral over time and the one 
over relative spatial positions provide a factor r^, since 
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(t) is different from zero only in a spacetime region 
of size r^, and the remaining integral over global spa- 
tial positions provides a factor the typical low-energy 
spatial volume. Putting everything together, we see that 
the diffusion term is of order with 

rji^ = 2ni(l + Si) — 2. This quantity defines the inverse 
of the decoherence time r^. Therefore, the ratio between 
the decoherence time and the low-energy length scale 
I is 

ra/l-e-^[Y.{h/ir-+^-\ (3.12) 

Only gravitational fluctuations whose size is very close 
to Planck length will give a sufficiently small decoher- 
ence time, because of the exponential dependence of 
e ~ e~'^^''-'/^(r//*)^. Slightly larger fluctuations will have 
a very small effect on the unitarity of the effective theory. 
For the interaction term that corresponds to the mass 
of a scalar field, the parameter rj vanishes and, conse- 
quently, Td/l ^ e^^. Thus, the scalar mass term will lose 
coherence faster than any other interaction. Indeed, for 
higher spins and/or powers of the field strength, rj > 1 
and therefore Td/l increases by powers of l/U^ For in- 
stance, the scalar-fermion interaction term c^^i/'V'j which 
has the next relevant decoherence time, corresponds to 
a decoherence ratio Td/l ^ e~'^l/U. We see that the de- 
coherence time for the mass of scalars is independent of 
the low-energy length scale and, for gravitational fluctu- 
ations of size close to Planck length, e may be not too 
small so that scalar masses may lose coherence fairly fast, 
maybe in a few times the typical evolution scale. Higher 
power and/or spin interactions will lose coherence much 
slower but for sufficiently high energies l^^, although 
much smaller than the gravitational fluctuations energy 
, the decoherence time may be small enough. This 
means that quantum fields will lose coherence faster for 
higher-energy regimes. 

IV. QUANTUM BATH 

In this section, we will take into account the quantum 
dynamical character of the fields a' that represent space- 
time gravitational fiuctuations (e.g., virtual black holes 
or quantum time machines) and describe spacetime foam 
in terms of a quantum thermal bath. By comparing the 
system consisting of low-energy fields suitably coupled to 
a quantum bath ||l^;|l^ with the results obtained above 
for gravitational fiuctuations, we will see that spacetime 
foam can be substituted by an effective quantum thermal 
bath. 

Let us start studying a system with a Hamiltonian 

H = Ho + Hint + Hi,, (4.1) 

where Hq is the bare Hamiltonian that represents the 
low-energy fields and _ffb is the Hamiltonian of a bath 



that, for simplicity, will be represented by a real massless 
scalar field. The interaction Hamiltonian will be chosen 
to have the form Hint — Chi, the noise operators ^ being 
given by 

£,^[x,t) = I f ^x-(^)[a+(fc)e*^"*"''"^ -a(fc)e-*("*-'=")]. 
J 

(4.2) 

In this expression, a and a+ are, respectively, the anni- 
hilation and creation operators associated with the bath, 
Lu = Vk^, and x-(i-li) are real functions that represent the 
coupling between the system and the bath for each fre- 
quency LU and for each interaction hi. These couplings 
X-(u;) can also be written in the position representation 
if we note that the momentum of the bath scalar field 
p{x, t) has the form 

p{x, t)=i J rffcV^[a+(fc)e*("'*-*=^) - a(A:)e-'('^*-'=^)], 

(4.3) 

so that the noise operators ^ have the form 

Cix,t) = jdx'x^{x-x')p{x',t). (4.4) 

Here, 

/dk 
-^X-{^)cos{ky) (4.5) 

represent the couplings between the low-energy field and 
the bath in the position representation. Since we are 
trying to construct a model for spacetime foam, we will 
assume that the couplings xKu) will be concentrated on 
a region of radius r and therefore the couplings x-(i^) 
will induce a significant interaction with all the bath fre- 
quencies u) up to the natural cutoff . Furthermore, 
these couplings have dimensions of length and we will 
also assume that they are of order e~^^'^'^^'^r. All the 
relevant information about the couplings is encoded in 
the commutation relations and the correlation function 
of the noise operators 

Let us start with the commutation relations at different 
times of the noise variables. Taking into account the 
commutation relations for the annihilation and creation 
operators a and , i.e., 

[a{k), a{k')] = [a+{k), a+{k')] = 0, (4.6) 
[a(fc),a+(fc')] = (5(fc- fc'), (4.7) 

it is easy to see that 

[C{t),eit')]^^fHt-t'), (4.8) 

where 
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/y(T)= / dujG'^{Lu)cos{LUT), (4.9) 

JO 

sm(uj\xi — xA) 
(c) = 8^^^ ^x-ico)xH^)- (4.10) 

Uj\Xi — Xj I 

Note that the functions G^^ (uj) and, hence, f'^-' {t) depend 
on the relative spatial distance \xi — Xj\, are symmetric 
in the pair of indices ij and are uniquely determined by 
the couplings X'i'-^) vice versa. In particular, they 
are completely independent of the state of the bath or 
the system. 

In order to compare this model with that of topologi- 
cal fluctuations previously described, it is convenient to 
introduce the so-called commutative noise representation 
jl3| by defining new noise operators a* in the following 
form: 

a\t)Q{t')^^[e{t),Q{t')]+ (4.11) 

for any operator Q. As we have seen, the commutators 
of the noise operators i^' at different times are c-numbers. 
Therefore, it is straightforward to check that the opera- 
tors a* commute at any time, i.e., 

[a\t),a^it')]=0. (4.12) 

However, the commutator of a' with any low-energy op- 
erator A is in general nonvanishing and has the form: 

[A{t),a'{t')]^ fdT[A{t),h,{T)]r{t' -t), (4.13) 

with hi{t) = U+{t)h,U(t) and U{t) = e"'^*. The func- 
tion f^^ (t) can be interpreted as a kind of memory func- 
tion. Indeed, these commutators are nonzero for low- 
energy operators that are in the future or, at most, in 
the near past of the noise and vanish only when they are 
in the far past. Only in the so-called first Markov ap- 
proximation the frontier among both regimes is sharply 
located where both noise and low-energy fields are at the 
same instant of time. 

We are now ready, following similar steps to those out- 
lined in the previous section, to write down the master 
equation for the low-energy density matrix. We will de- 
scribe the whole system (low-energy field and bath) by a 
density matrix pT{t). We will assume that, initially, the 
low energy fields and the bath are independent, i.e., that 
at the time to 

p^ih) ^ p{ta) ® Ph- (4.14) 

As in the classical noise case, if the low-energy fields and 
the bath do not decouple at any time, an extra renormal- 
ization term should be added to the Hamiltonian. In the 
interaction picture, the density matrix has the form 

p\{t) = U+{t)pr{t)U{t), (4.15) 



with Uit) = C/o(i)C/b(i), where Ua{t) = e-'""* and 
Uh{t) = e^'^t,*. It obeys the equation of motion 

p\it) = -i[em{t),p\it)]. (4.16) 

Here, 

at) = u+{t)cu{t) = u+{t)eu^{t), (4.17) 

h\{t) = U+{t)h,U{t) = U+{t)h,Uo{t). (4.18) 

Integrating this evolution equation and introducing the 
result back into it, we obtain the following integro- 
differential equation: 

p\{t)^-z[eimit),p\{to)] 

- fdt'[e{t)h\{t)A^'{t')h){t').p\{t')]]. (4.19) 

Jto 

If we now trace over the variables of the bath, define 
p\t) EE trb[p!r(t)] and note that tr^l^' {t)h\{t) p\{to)] = 
(because trb[C*(t)pb] = 0), we obtain 

p\t) = - fdt'tr^ {[CmMt), [e{t')h},{t').p\{t')]]) . 

Jto 

(4.20) 

In the weak-coupling approximation, which implies 
that ^'/li is much smaller than Hq and 7?b (this is jus- 
tified since it is of order e), we assume that the bath 
density matrix does not change because of the interac- 
tion, so that p\{t) — p^{t) (X) pb- The error introduced by 
this substitution is of order e and ignoring it in the mas- 
ter equation amounts to keep terms only up to second 
order in this parameter. Since h^j{t')] — because 

[Cihj] — 0, the right hand side of this equation can be 
written in the following way 

- 1 fdt'{{[C{t),i'{t')]+)[h\{t), [h){t'),p\t')]] 

H[c{t).e{n)[h\{t)Ah){t'),p\t')u]), (4.21) 

where the average of any operator Q has been defined 
as (Q) = trb((3pb)- Next we note that {[Cit),i^ {t')]) = 
ip^{t — t') and, using the commutative noise representa- 
tion introduced above, we can write 

^([f W,e(0]+> = {a\t)a^t'))^c^^{t-t'). (4.22) 

If we make the assumption that the bath is in a ther- 
mal state ph = Z~^e~^^/'^ with a temperature inversely 
proportional to the size of the gravitational fluctuations 
(e.g., the radius of the virtual black holes, or the size of 
the regions containing closed timelike curves in the case 
of quantum time machines), T ^ l/r, the correlation 
function {t — t') acquires the form: 

/•oo 

c*J'(r)=/ dwijjG'3{uj)[N{ijj) + l/2]cos{ujT), (4.23) 
Jo 
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where N{oj) = [exp(w/T) — is the mean occupation 
number of the bath corresponding to the frequency u>. 
In this calculation, we have made use of the following 
relations, valid for a thermal state: 

{a{k)) = {a+{k))=0, (4.24) 
(a(fc)a(fc')) = {a+{k)a+{k')) = 0, (4.25) 
{a+{k)a{k')) = N{oj)5{k - k'). (4.26) 

Similarly, we can easily compute the higher order corre- 
lations {a^ {t)a^ {t')a'' it")) , etc. Those containing an odd 
number of fields a' turn out to be identically zero while 
those containing an even number can be written in terms 
of the two-point correlation function (t) . This means 
that the trace {Q) corresponds to a Gaussian average 
over a*, provided that the bath is in a thermal state, as 
we are considering. In this way, we have established a 
relation between a quantum thermal bath and spacetime 
foam, which can also be described by a Gaussian average, 
as we have seen. 

The Markov approximation allows the substitution of 
p'(t') by p'(i) in the master equation because the inte- 
gral over t' will got a significant contribution from times 
t' that arc close to t due to the factors f^^it — t') and 
d^it — t') and because, in this interval of time, the den- 
sity matrix will not change significantly. Indeed, the 
typical evolution time of is the low-energy time scale 
/, which will be much larger than the time scale r as- 
sociated with the bath. If we perform a change of the 
integration variable from t' to t = t — t', write 

p\t') = p\t -r)= p\t) - Tp\t) + 0{r^), (4.27) 

and introduce this expression in the master equation 

above, we easily see that the error introduced by the 
Markovian approximation is of order e^, i.e., it amounts 
ignore a term of order e**. The upper integration limit 
t in both integrals can be substituted by oo for evolu- 
tion times t — to much larger than the correlation time r, 
because of the factors (r) and c*^ (r) that vanish for 
T > r, which is equivalent to take the initial condition to 
the infinite past to — > — oo. 

Then, the master equation in the interaction picture 
acquires the form 

i r°° ■■■ 

^'^^^ = -2 I drf'i^)WM[h){t-T),p\t)]+] 

/•OO 

- / drd\T)m),[h)it-r),p\m. (4.28) 
Jo 

We can now transform the resulting equation back to the 
Schrodinger picture 

i r°° ■ 

p = -i[Ho,p] - - rfrr (r)[/i,, [h){-T),p]+] 

- / dTc'^{r)[hi,[h){-T),p]]. (4.29) 
Jo 



After an integration by parts, the second term of the 
right hand side becomes 

/•oo 

'-P^{0)[h,h„p]-- dTr{T)[h,,[h){-T),p]+]. 

(4.30) 

The first term is just a finite renormalization of the orig- 
inal low-energy Hamiltonian from Hq to 

Hi, = Ho - ^fHO)hihj (4.31) 

and the master equation can then be written in its final 
form 

p=-t[Hip]~- dTr{T)[K,m-T),p] + ] 

POO 

- / dTc'^{T)[hi,[h){-T),p]]. (4.32) 
Jo 

Before discussing this equation in full detail, let us first 
study the classical noise limit. With this aim, let us 
introduce the parameter 

a = j dk'[a{k),a+{k% (4.33) 

which is equal to 1 for quantum noise and for clas- 
sical noise. Then, the /-term is proportional to a and 
therefore vanishes in the classical noise limit. The c-term 
also contains a factor a but, in addition, N{(jj) becomes 
N(cruj) when introducing the parameter a. In the limit 
cr ^ 0, the term proportional to 1/2 in c*-^ (t) vanishes 
and the term proportional to N{aco) acquires the value 
c^i^^sir) = TP^{t). Also, the renormalization term of 
the low-energy Hamiltonian vanishes in this limit. In 
this way, we have arrived at the same master equation 
that we obtained in the previous section. This is not sur- 
prising because the origin of the /-term is precisely the 
noncommutativity of the noise operators, i.e., its quan- 
tum nature, while the Cdass-term actually contains the 
temperature effects. At zeroth order in r/l, the master 
equation for classical noise then acquires the form 

p = -i[Ho, p] - / dTcl,,(r) [hi, [hj,p\\ . (4.34) 
Jo 

Let us now analyze the general master equation, valid 
up to second order in e that takes into account the 
quantum nature of the gravitational fluctuations. These 
contributions, although small in the low-energy regime, 
might still be experimentally testable. In addition, they 
may provide interesting information about the higher- 
energy regimes in which I may be of the order of a few 
Planck lengths and for which the weak-coupling approx- 
imation is still valid. In order to see these contributions 
explicitly, let us further elaborate the master equation. 
In terms of the operator Lq defined as LqA = [Ho,A\ 
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acting of any low-energy operator A, the time dependent 
interaction r) can be written as 

h){-T) ^e~'^°''hj. (4.35) 

The interaction hj can be expanded in eigenoperators 
of the operator Lo, i.e., 

h, - J d^in (/i+o + hjn) , (4.36) 

with Lf)h^^ — ztflh^Q and d^^ being an appropriate 
spectral measure, which is naturally cut off around the 
low-energy scale This expansion always exists pro- 
vided that the eigenstates of Hq form a complete set. 
Then, h^j{—T) can be written as 

h)i-r) - J dt,n{e-''''h+^ + e'^^/jT^). (4.37) 

It is also convenient to define the new interaction opera- 
tors for each low-energy frequency f2 

h]n - h+, - h-^, (4.38) 
h'n = h+n + hjn- (4.39) 

Both the term proportional to /'-'(r) and the term 
proportional to d^r) are integrated over r £ (0,cx)). 
Because of these incomplete integrals, each term provides 
two different kinds of contributions coming from the bulk 
term and the principal part in the well-know formula 

/>OC 

/ dre"^^ ^TTS{uj)+V{i/uj), (4.40) 
Jo 

where V is the Cauchy principal part . 

The master equation can then be written in the follow- 
ing form 

P = -(i^O + -^diss + LdiS + iLstark + «ilamb)p, (4.41) 

where the meaning of the different terms are explained 
in what follows. 

The first term —iL'^p^ with L'^p — [Hq^ p], is re- 
sponsible for the renormalized low-energy Hamiltonian 
evolution. The renormalization term is of order as 
compared with the low-energy Hamiltonian Hq, where 

= Y.ij{^*l^f-'^'^- and, remember, 77, = 2niil+Si)-2 
is a parameter specific to each kind of interaction term 
hi. 

The dissipation term 

Ldi..p = d^in^G''mh^, [h]n,p]+] (4.42) 

is necessary for the preservation in time of the low-energy 
commutators in the presence of quantum noise. As we 
have seen, it is proportional to the commutator between 
the noise creation and annihilation operators associated 



with the effective bath that represents spacetime foam 
and, therefore, vanishes in the classical noise limit. Its 
size is of order e^r/P. 

The diffusion process is governed by 

(4.43) 

which contains two contributions: the first one is a tem- 
perature effect of order /I and the second is a cold dif- 
fusion originated in the vacuum fluctuations of the grav- 
itational field and it is of order e^r/l'^. In the classical 
noise limit, only the first contribution survives and was 
already studied in the previous section. 

The next term provides an energy shift which can be 
interpreted as a gravitational ac Stark effect by compar- 
ison with its quantum optical analog ||l^,|l^ . Its expres- 
sion is 

//>oo /-J 
dpnVj^ diu-^^G'^u;)N{co)[h,, [h^n.p]]- 

(4.44) 

Although it is also a temperature-dependent effect with 
the same origin as the diffusion term, it contains a 
Cauchy principal part. This translates into the fact that 
it is smaller than the diffusion term although it does not 
vanish in the classical noise limit. It is of order e^r/P. 

Finally, LiambP is an energy shift generated by the vac- 
uum fluctuations of the gravitational field (as the dissi- 
pation term and the cold diffusion term) and that can 
therefore be interpreted as a gravitational Lamb shift. It 
has the form 

If r°° n 

LiamhP = dpnV doj^—^C-^iuj) 

X {,o[h,, [h],„p]] - n[h,, [h'jn,p]+]} ■ (4.45) 

The second term is of order e'^r'^ /l^, which is fairly small. 
However, the first term will provide a significant contri- 
bution of order e'^r/l'^\og{l/r). This logarithmic depen- 
dence on the relative scale is indeed characteristic of the 
Lamb shift ||,|l|,||. 

As a summary, the c-term gives rise to four different 
contributions: a thermal diffusion term, another diffu- 
sion term originated from the vacuum fluctuations of the 
bath, a contribution to what can be interpreted as a grav- 
itational Lamb shift, and, finally, a shift in the scalar- 
field oscillation frequencies that can be interpreted as a 
gravitational Stark effect. The /-term provides a dissipa- 
tion part, necessary for the preservation of commutators, 
and another contribution to the gravitational Lamb shift. 
The size of these effects generated by spacetime foam, 
compared with the bare evolution, are the following: the 
thermal diffusion term is of order e^, which is the only 
one that survived in the approximations of the previous 
section; the diffusion created by vacuum fluctuations, the 
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damping term, and the Stark effect are smalier by a fac- 
tor r/l] and the Lamb shift has two contributions: one 
is smaller than the diffusion term by a factor {r/lY and 
the other is of order {r/l) \og{l/r) as compared with the 
diffusion term. Note that the quantum effects induced by 
spacetimc foam become relevant as the low-energy length 
scale I decreases, as we see from the fact that these ef- 
fects depend on the ratio r/Z, while, in this situation, the 
diffusion process becomes slower, except for the mass of 
scalars, which always decoheres in a time scale which is 
close to the low-energy evolution time. 

These quantum gravitational effects could be mea- 
sured, at least in principle, since they are just energy 
shifts and decoherence effects similar to those appearing 
in other areas of physics, where fairly well established ex- 
perimental procedures and results exist, and which can 
indeed be applied here — such as those briefly discussed 
below — provided that sufficiently high accuracy can be 
achieved. On the other hand, scalar fields lose quan- 
tum coherence extremely fast and Hawking has argued 
H] that this might be the reason for not observing the 
Higgs particle. He has also suggested that loss of quan- 
tum coherence might be responsible for the vanishing of 
the angle in quantum chromodynamics 

Neutral kaon beams have been proposed as experi- 
mental systems for measuring loss of coherenc e owing to 
quantum gravitational fluctuations 
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In these 

systems, the main experimental consequence of the diffu- 
sion term (together with the dissipative one necessary for 
reaching a stationary regime) is violation of CPT ||2^ , p^ 
because of the nonlocal origin of the effective interac- 
tions. The estimates for this violation are very close to 
the values accessible by current experiments with neu- 
tral kaons and will be within the range of near-future 
experiments. Macroscopic neutron interferometry ||r^ , |27|| 
provides another kind of experimental systems in which 
the effects of the diffusion term have measurable conse- 
quences since they may cause the disappearance of the 
interference fringes |l5| , |27| . 

As for the gravitational Lamb and Stark effects, they 
are energy shifts that depend on the frequency, so 
that different low-energy modes will undergo different 
shifts. This translates into a modification of the disper- 
sion relations, which makes the velocity of propagation 
frequency-dependent, as if low-energy fields propagated 
in a "medium". Therefore, upon arrival at the detec- 
tor, low-energy modes will experience different time de- 
lays (depending on their frequency) as compared to what 
could be expected in the absence of quantum gravita- 
tional fluctuations. These time delays in the detected 
signals will be very small in general. However, it is still 
possible to measure them if we make the low-energy par- 
ticles travel large (cosmological) distances. In fact, 7- 
ray bursts provide such a situation as has been recently 
pointed out |^ , thus opening a new doorway to observa- 
tions of these quantum gravitational effects. Indeed, the 
ratio between the time delay owing to gravitational fluc- 
tuations and the width of the intrinsic time structure of 



7-ray bursts has been estimated to be of order 1 for emis- 
sions with millisecond time structure and energy around 
20 MeV , provided that they travel a distance of 10^" light 
years |16| , which are compatible with 7-ray burst obser- 
vations. If this sensitivity can actually be reached, one 
would expect that the presence of the gravitational Lamb 
and Stark shifts predicted above could be observationally 
tested. 



V. VIRTUAL BLACK HOLES, WORMHOLES, 
AND TIME MACHINES 

It is well-known that it is not possible to classify all 
four-dimensional topologies Q and, consequently, all the 
possible components of spacetime foam. Here, we will 
briefly discuss three different kinds of fluctuations: sim- 
ply connected nontrivial topologies, multiply connected 
topologies with trivial second homology group (i.e. with 
vanishing second Betti number), and finally spacetimes 
with a nontrivial causal structure, i.e., with closed time- 
like curves, in a bounded region. 

The effective description proposed in this paper and 
the associated master equation are particularly suited to 
the study of low-energy effects produced by simply con- 
nected topology fluctuations (e.g., virtual black holes). 
Hawking |6| has shown that compact simply connected 
bubbles with the topology S"^ x S'^ (whose second Betti 
number is i?2 = 1) can be interpreted as closed loops of 
virtual black holes if one realizes that the process of 
creation of a pair of real charged black holes accelerating 
away from each other in a spacetime which is asymptotic 
to 3?^^ is provided by the Ernst solution This solu- 
tion has the topology x 5^ minus a point (which is 
sent to infinity) and this topology is the topological sum 
of the bubble S"^ x S"^ plus SR^. Virtual black holes will 
not obey classical equations of motion but will appear as 
quantum fiuctuations of spacetime and thus will become 
part of the spacetime foam. Particles could fall into these 
black holes and be re-emitted. The scattering amplitudes 
of these processes |^ could be interpreted as being pro- 
duced by nonlocal effective interactions that would take 
place inside the fluctuations and the master equation ob- 
tained above could then be interpreted as providing the 
evolution of the low-energy density matrix in the presence 
of a bath of ubiquitous quantum topological fiuctuations 
of the virtual-black-hole type. 

Wormholes Q, i.e., multiply connected fiuctuations 
(with vanishing second Betti number), also admit a de- 
scription in terms of nonlocal interactions that, in the 
weak-coupling approximation, become bilocal. These 
quantum fiuctuations connect spacetime points that may 
be far apart from each other, in the dilute gas approx- 
imation. Therefore, diffeomorphism invariance on each 
spacetime region requires the coefficients of this bilo- 
cal interaction term to be spacetime independent. The 
same conclusion can also be reached if we analyze worm- 
holes from the point of view of the universal covering 
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manifold, which is, by definition, simply connected. A 
wormhole is then represented in the universal covering 
manifold by two boundaries, suitably identified, located 
at infinity. This identification can be implemented by 
introducing coefficients c*-' that relate the bases of the 
Hilbert space of wormholes in both regions of the univer- 
sal covering manifold. As coefficients in a change of basis, 
c*-' cannot depend on spacetime positions and, therefore, 
will just be constant. This means that the fields a* can- 
not be interpreted as noise sources that are Gaussian dis- 
tributed at each spacetime point independently, because 
the correlation time for the fields a' is infinite. Indeed, 
the constancy of implies that they are infinitely coher- 
ent and the Gaussian distribution to which they arc sub- 
ject is therefore global, spacetime independent l^. One 
could still expect some effects originated in the quan- 
tum nature of a* such as a cold diffusion term in the 
master equation or even dissipation. However, because 
they are spacetime independent, they commute with ev- 
ery operator, including low-energy ones, thus giving rise 
to superselection sectors. Therefore, all the terms in the 
master equation, except the one responsible for the uni- 
tary low-energy evolution, vanish. Still, wormholes can 
be represented by a thermal bath as we have done with 
localized gravitational fluctuations. However, in order to 
reproduce their infinite correlation time, the couplings 
^ between the bath and the low-energy fields must be 
constant, they must commute with every other operator, 
and, related to these two facts, only the zero-frequency 
(i.e., infinite wavelength) mode of the bath can be cou- 
pled to the low-energy fields, thus leading to a unitary 
effective theory. 

From the semiclassical point of view, most of the hith- 
erto proposed time machines Q are unstable because 
quantum vacuum fluctuations generate divergences in the 
stress-energy tensor, i.e. are subject to the chronology 
protection conjecture [ pl[ . However, quantum time ma- 
chines confined to small spacetime regions, for which 
the chronology protection conjecture does not apply 
are likely to occur within the realm of spacetime foam, 
where strong causality violations or even the absence of a 
causal structure are expected. These Planck-size regions 
with quantum time machines admit an effective represen- 
tation in terms of nonlocal interactions that account for 
the causality violations and will lead to a loss of quan- 
tum coherence [Q that can also be effectively described 
as coming from the interaction of low-energy fields with 
a thermal bath. In this case, the low-energy density ma- 
trix will also evolve according to the master equation 
obtained in the previous sections. 



VI. CONCLUSIONS 

In this paper, we have built an effective theory in 
which quantum gravitational spacetime foam has been 
substituted by a fixed classical background plus nonlo- 



cal interactions between the low-energy fields confined 
to bounded spacetime regions of nearly Planck size. In 
the weak-coupling approximation, these nonlocal inter- 
actions become bilocal. The low-energy evolution is not 
unitary because of the absence of a nonvanishing timelike 
Hamiltonian vector field. The nonunitarity of the bilo- 
cal interaction can be encoded in a classical noise source 
locally coupled to the low-energy fields and subject to 
a Gaussian probability distribution. Then, the evolu- 
tion of low-energy fields is provided by a master equa- 
tion which contains a diffusion term. This diffusion is a 
direct consequence of the nonlocal character of the quan- 
tum gravitational fluctuations encompassed by spacetime 
foam. The decoherence rate is suppressed by powers of 
the ratio between the gravitational fluctuation size and 
the low-energy length scale, except for the mass interac- 
tion term of scalar fields for which this rate is comparable 
with the low-energy evolution scale. 

We have argued that the quantum nature of spacetime 
foam is not represented in this effective theory but only 
its thermal properties. A model in terms of a quantum 
thermal field, which in the classical noise limit coincides 
with the one described above, has been proposed as de- 
scribing the quantum and thermal properties of space- 
time foam. In this model, the low-energy density matrix 
evolves according to a master equation that, apart from 
inducing loss of coherence, contains additional terms that 
may be relevant for sufficiently high energies. These 
terms correspond to a dissipation process that ensure the 
preservation of commutators, a cold diffusion, and energy 
shifts that can be interpreted as gravitational Lamb and 
Stark effects. We have also briefly discussed some of the 
possible experimental implications that these quantum 
gravitational effects may have. A constructive model in 
terms of nonlocal interactions that takes into account 
the quantum origin of spacetime foam will be developed 
elsewher e |3^j3^ within the formalism of Feynman and 
Vernon iJ^. 

Finally, among the possible components of spacetime 
foam, the role of virtual black holes and small bounded 
regions that contain closed timelike curves has been 
briefly analyzed in the context of our effective model. 
We have also argued that dilute wormholes do not admit 
a description in terms of a thermal bath coupled to the 
whole low-energy spectrum and that they are inflnitely 
coherent as was already shown by Coleman pi . 
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